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PREFACE 


The work described in this publication was sponsored 
by the United States Army Intelligence Center and School. 
The writing and publication of this paper was supported 
by the Jet Propulsion Laboratory, California Institute 
of Technology, under a contract with the National 
Aeronautics and Space Administration, NAS 7-918, RE 182, 
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EXECUTIVE SUMMARY 


This Technical Memorandum was prepared 
originally as part of the Generic ELINT/COMINT Sensor 
Report (FY-86) which was elmininated under the FY-87 
statement of work (SOW #2) , undated (delivered to 
JPL 19 November 1986) . 

The purpose of the Generic ELINT/COMINT Sensor 
Report, of which this paper was intended (in its 
final form) to become part of, was to establish a 
basic superhetrodyne receiver based sensor model and 
perform simulations with it to determine the shaping 
or coloring of the statistical distributions of the 
radar free-space signal parametrics by a typical 
sensor prior to reaching the self-correlation 
processes. It was also intended for incorporation 
into the algorithm test bed so algorithms could be 


tested with realistic 
unrealistic stastically 

distorted 
pure data . 

data 

rather 

than 

This 

work was 

originated 

in 

support 

of 


unanswered questions from previous self-correlation 
studies. The modeling and simulation approach was 
used because "live date" could not be obtained. 

This paper is being published because it was 
completed in FY-86 with FY-86 funds and still serves 
a useful function. 

Several significant results were noted while 
performing simulations on this model. The first 
order math model is developed for the Bragg-cell 
receiver. Whether the input distributions were 
uniform or Gaussian, the output distributions were 
found to be shaped so as to append a staircased 
triangular distribution. 
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ABSTRACT 


A 

Bragg-cell receivers are employed in specialized Electronic Warfare 
(EW) applications for the measurement of frequency. Bragg-cell receiver 
characteristics are fully characterized for simple RF emitter conditions, 
but less understood for complex and wideband RF emitter signals. This 
receiver is early in its development cycle when compared to the IFM receiver. 

Functional mathematical models are derived and presented in this 
report for the Bragg-cell receiver. Theoretical analysis is presented and 
digital computer signal processing results are presented for the Bragg-cell 
receiver. Probability density function analysis are performed for output 
frequency. 

Probability density function distributions are observed to depart from 
assumed distributions for wideband and complex RF signals. This analysis 
is significant for high resolution and fine grain EW Bragg-cell receiver 


systems. 

■A 



NEW TECHNOLOGY 


None 
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INTRODUCTION 


The development effort is for theoretical analysis, software algorithms 
implementations, signal processing analysis, and test analysis for probability 
density function characterization of a Bragg-cell receiver system. SigPro 
Systems Inc. is performing this work under JPL Contract No. 957474, Mod. 
No. 1. 

This final report describes the completed study on the Bragg-cell 
receiver. A functional mathematical model is derived for the receiver. 
This math model is implemented into digital computer software and signal 
processing analyses are performed. 

Theoretical analysis is developed for the Bragg-cell receiver. Ideal 
and practical subsystem implementations are presented and analyzed. The 
key ELINT parameter of interest is frequency, hence, the theoretical analysis 
presents the basic analysis approach for frequency. Probability density 
function analysis is presented for first-order subsystem model. 

The theoretical model for the Bragg-cell is developed into a digital 
computer program. All computer programs are listed in Appendix A. 

Signal processing analysis is performed using the derived first-order 
model. The pdf characteristics of the Bragg-cell receiver’s output frequency 
are determined for assumed input frequency pdf (constant and Gaussian) 
characteristics. 
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BRAGG CELL RECEIVER 


INTRODUCTION 

A Bragg-cell receiver is used to measure frequencies of input RF radar 
pulses, with wide frequency coverage, large dynamic ranges and simultaneous 
or time coincident signals. The ideal Bragg-cell receiver will measure the 
frequency parameter on a single pulse basis with no distortion. Practical 
Bragg-cell receivers will distort the frequency characteristic of a radar emitter 
because the ideal Bragg-cell receiver characteristics are difficult to approximate 
in real life. 

The radar's frequency parameter can be characterized using a random 
variable. The random variable may be quantified by the probability density 
function (pdf), mean value, variance and other moments. As a radar signal 
passes through the Bragg-cell receiver the random signal description is modified 
or distorted by the receiver's nonlinear transfer function characteristic. 
Bragg-cell receiver characterization of RF frequency will differ from actual 
radar emitter characterization by multifaceted distortion effects. 

This section presents a functional description of the Bragg-cell receiver 
with potential distortion sources and mechanisms and a first-order Bragg-cell 
receiver model is developed. This Bragg-cell receiver model is used to analyze 
the distortion of the frequency's pdf. The first-order Bragg-cell receiver model 
is limited to key distortion effects. 

IDEAL BRAGG-CELL RECEIVER CHARACTERIZATION 

A functional block diagram of a Bragg-cell receiver is shown in Figure 
1. The basic receiver includes an RF-to-IF downconverter, Bragg-cell with 
AO transducer, laser source, beam expander optics, beam focusing optics, 
detector array and processing electronics. 
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The RF input signal is characterized by 


— XC'fc) cosT^-t + 


E q. 1 


where 

x(t) = Pulsed or CW signal, 

amplitude modulation function, 

UJ C = fixed angular carrier frequency, 

= phase deviation function. 

The RF-to-IF downconverter is described in Reference 1. The ideal 
downconverter is composed of an ideal local oscillator which is a perfect CW 
signal, a perfect product mixer, and an ideal bandpass filter. The ideal IF 
signal output from the downconverter is described quantitatively by: 


ijCt") = XCt) cos^ujiyt + 4’ft'T] 


Eq. 2 


Equation 2 is identical to Equation 1 with the fixed angular frequency being 
changed from to w I5 ,. The IF center frequency is chosen to match the 
center frequency of the Bragg-cell and AO transducer. Ideally the IF signal 
replicates the input RF signal. The RF-to-IF downconverter is described and 
characterized in Reference 1. RF-to-IF downconverter analysis will not be 
presented in this report. 

The IF signal drives the AO transducer which launches an acoustic signal 
through the Bragg-cell. Specifically, the IF electronic signal is converted 
into an acoustic wave (slow wave replication of the IF signal) which propagates 
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through the optically transparent Bragg-celL, Through the elasto-optic effect, 
the acoustic wave produces a spatial modulation of the refractive index in 
the Bragg-cell (Reference 2.) 

The laser source is expanded by the beam expander optics and illuminates 
the Bragg-cell. As the laser light passes through the Bragg-cell, the refractive 
index variations produced by the acoustics wave (which is a slow form replica 
of the IF signal) are impressed onto the optical signal as a spatial phase 
modulation. 

Focusing optics are used to equivalently Fourier transform the modulated 
optical beam. Resultant optical Fourier transform signal is focused onto the 
detector array or frequency focal plane (see Figure 1.) The photodiode detector 
array is used to capture and convert the optical signal to an electronic signal 
representation. The detector array is a linear photodiode array which is 
functionally used to detect and measure light intensity versus frequency of 
the transformed optical signal. The detector array is electronically read by 
the processing electronics, which also provides follow-on signal processing 
functions. 

The Bragg-cell modulator is depicted in Figure 2. The incident laser 
light interacts with the acoustic signal in the cell to produce a diffracted light 
information signal along with undeflected laser light signal. The diffracted 
laser signal is also frequency shifted by the acoustic signal frequency. The 
incident laser light hits the Bragg-cell at an angle 9^ with respect to the 
a axis. The IF electronic signal is sent to the AO transducer to convert the 
IF signal to a replica acoustic signaL The acoustic wave of frequency f 1 
propagates with velocity along the x - axis in the Bragg-cell. 

The diffracted laser light signal results from the AO interaction effects 
between the incident laser beam and the acoustic information bearing signal 
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in the Bragg-cell. Maximum diffracted laser light energy occurs in the first 
order beam when (Reference 2): 



sin' 


e vT 


and 



Eq. 3. 


Eq. 4 


Wavelength of laser light 
Frequency of acoustic wave 
Velocity of acoustic wave in Bragg-cell along the 
x-axis. 

Effective angle of incident laser light. 

Effective angle of first order diffracted laser light 
information signal. 

Detailed Bragg-cell analysis would require AO cell refractive index 
variations analysis for precise quantification. This first order Bragg-cell model 
does not require this detailed analysis. Functionally the Bragg-cell model 
is accurate and conceptually clear for the external angles 9 j and 9^. 

Bragg-cell processors are, in general, configured so that 9.« 0.1 
radian, hence. Equation 3 can be readily and accurately simplified to: 


where 


Vs 


© •» - 
a, = 


10 



Eq. 5. 



Using Equation 4 and 5, the foEowing equation for the effective angle 
of the first order diffracted laser light information signal is: 



ZJrx 

2 v s 


Eq. 5 


Equation 6 revesils that 9 d varies directly with the acoustic wave frequency. 
As stated earlier, the IF signal frequency is related to the acoustic wave 
frequency by: (Reference 2) 



Eq. 7. 


where K is a constant. 

Equation 8 results by combining Equations S and 7, 




Eq. 8. 


Equation 8 shows that 9 d varies directly with the IF signal carrier frequency. 

The actual displacement of the diffracted laser information signal from 
the undeflected laser light (see Figure 2) is given by (Reference 2.) 



F 3- 

K Vg 


Eq. 9. 
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where F is the focal length of the focusing optics. 

A fixed Bragg-cell receiver will result in Equation 9 reducing to: 

d — K \ ‘^'XF Et? * 10 

The actual displacement, d, of the diffracted laser beam relative to 
the undeflected laser light passing through the Bragg-cell is directly proportional 
to the IF frequency of the input IF electronic signal. The diffracted laser 
beam will hit the linear detector array at a specific location which corresponds 
to one IF frequency. The Bragg-cell receiver is calibrated by varying the IF 
frequency of the input IF signal over the entire operating bandwidth and 
capturing and analyzing the electrical output signal from the linear detector 
array. Detector array element signal versus IF frequency characteristic is 
used to quantify the Bragg-cell receiver's frequency measurement performance. 

The undeflected laser light is (the zero-order output from the Bragg-cell) 
described by: 

JL 0 ($ = Ufe) Cas 

where 

-&«(£)* undeflected laser light signal. 

L,($ * amplitude or intensity of undeflected laser light signal. 

Tji = laser light frequency 

OrCf) = phase variations on laser light signal. 

The diffracted laser information signal is quantified mathematically 
as: 


+■ OrC-fcfJ 


Eq. 11. 
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Eq. 12. 



L/t) Cos 



where 

Diffracted laser information signal. 

Lj(t') = Amplitude or intensity of diffracted laser information signal. 

? t = Frequency of information signal. 

%Lt) = Phase variations of information signal. 

Equation 12 reveals that the diffracted laser output signal is frequency 
shifted. The frequency shift is produced by the laser light signal and acoustic 
signal interaction in the Bragg-celL 

The Bragg-cell receiver is conceptually a parallel, multi-channel spectrum 
analyzer, which determines the Fourier magnitude spectrum of the input IF 
signal. The electronic circuit model for the Bragg-cell receiver is depicted 
in Figure 3. The IF signal input is applied to a parallel narrowband filter bank 
with an associated diode detector and filter bank. Each output signal is a 
power or energy indicator of the IF signal spectral content in each selected 
narrowband filter. 

PRACTICAL BRAGG-CELL RECEIVERS 
The ideal Bragg-cell receiver is characterized in the previous section, 
with the measured frequency being derived from Equations 9, 10, and 12. 
This section will identify and briefly present some key differences between 
practical Bragg-cell receivers and the ideal bragg-cell receiver. The Bragg-cell 
receiver, shown in Figure 1, is a reference for this discussion. 
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Figure 3. 


Electronic Circuit Model for Bragg-cell Receiver 



BRAGG CELL 


The ideal Bragg-cell is assumed to have an undistorted acoustic signal 
replica (at acoustic frequency) of the IF input signal propagating through the 
cell. The AO transducer provides this electronic IF signal to acoustic signal 
conversion. The AO transducer has transfer function nonlinearities, which 
produces unwanted and added acoustic signal components. Simple sinusoidal 
electronic signals are subject to minor variations, while complex and wideband 
signals may be significantly distorted by the nonlinear transducer action. 

The ideal Bragg-cell is assumed to provide a perfect acoustic signal 
termination, at the far end of the cell. In general, the acoustic signal 
termination is not perfect, which cause low level acoustic signal reflections 
and distortion generated signal components to propagate back through the 
cell. These reflected acoustic signals are considered as additional distortion 
sources, hence, variations in the diffracted laser information signal can be 
expected. 

l 

Actual distortion effects are determined by AO transducer nonlinearities, 
reflected acoustic signal levels, total acoustic signal parameter descriptors 
and Bragg-cell properties. The complete acoustic signal's instantaneous 
frequency is determined by amplitudes, frequencies, and phases of desired 
acoustic signal, undesired reflected acoustic signals and nonlinearly generated 
and unwanted acoustic signals. Frequency variations produced by instantaneous 
frequency changes will produce angle variations in the diffracted laser 
information signal; thus, frequency measurement errors in the Bragg-cell 
receiver. 
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LASER SOURCE 


The laser light source is ideally assumed to be constant amplitude and 
fixed frequency signal. Actual laser sources can have some minor frequency 
shifts. Minor frequency shifts will results in small angle variations in the 
diffracted laser information signal. (See Equation 3.) 

Frequency variations will widen the effective spotsize and change the 
angle of the diffracted laser beam, which means the light energy may be spread 
over more detector elements in the linear detector array. 

The laser light signal does not have a constant intensity or amplitude. 
The amplitude is quantified by deterministic and random signal components. 
These amplitude variations will result in diffracted laser information signal 
intensity level variations at the detector plane. 

OPTICS 

Beam expander optics and diffraction signal focusing optics are 
fundamental modules in the Bragg-cell receiver. Beam expander optics shape 
the coherent laser beam to illuminate one entire side of the Bragg-cell. This 
beam expansion provides a long path interaction between the acoustic signal 
and laser light within the Bragg-celL Beam expander optics and laser source 
are positioned so that the laser light strikes the cell at the Bragg angle. 
Nonlinearities, positional variations, and rotational variations in the beam 
expander optics can produce amplitude and angle variations in the diffracted 
laser information signal. These errors are usually very small for well-designed 
optics. 

Focusing optics are used as the so-called Fourier transform lens. The 
Bragg-cell is positioned in the front focal plane of the lens. The photodetector 
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array is positioned at the back focal plane of the lense. (Reference 3.) Again, 
nonlinearities and location variations in the focusing optics variations can 
produce intensity and angle variations in the diffracted laser information signal. 
These errors are also very small in well-designed optics. 

DETECTOR ARRAY 

The ideal detector array is assumed to be a long linear array with high 
frequency resolution and sensitivity. Practical Bragg-cell receivers use long 
linear photodiode arrays, where typically 512 or 1024 photodiodes are closely 
spaced along the array length. 

Figure 4 presents an illustration of two diffracted laser information 
signals striking the photodetector array at two different diodes (two different 
frequencies.) In actual practice, the light spotsize for each beam is larger 
than the area of one photodiode; hence, the light spills onto the area between 
detectors and also onto adjacent detectors. 

Diffracted laser information signal angle variations (frequency variations) 
produce energy spread over several or many photodiode cells. Each photodiode 
cell is used to determine the signal energy in resolved and calibrated frequency 
space. The photodetector array is electronically readout to determine frequency 
and amplitude information for each diffracted laser light information signal. 
The photodetector array is functionally and practically described as a frequency 
sampling unit, with 512 to 1024 frequency bins. Each photodiode determines 
the elemental frequency resolution. 

Bragg-cell receivers can accurately measure frequencies of time coincident 
or time overlapped signals, which is a big advantage over other receivers, 
such as the EFM receiver. These receivers are available with bandwidths to 
approximately 1 GHz, and frequency resolutions of 100 KHz to 10 MHz. The 
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Figure *1. Simplified Diagram of Bragg-cell Receiver with Two Input Signals. 


Bragg-ceil receiver's dynamic range is approximately 30 to 40 db. '.vith the 
photodetector array's dynamic range being the critical limiting factor. 


PROBABILITY DENSITY FUNCTION ANALYSIS 
BRAGG-CELL RECEIVER 


The probability density function (pdf) of the Bragg-cell receiver's measured 
frequency parameter is considered in this section. The pdf characteristic 
is developed for a first order Bragg-cell receiver model under high 
signal-to-noise ratio conditions. The IF signal's frequency characteristic is 
assumed to have a constant pdf or a Gaussian pdf. 

Ideal and practical Bragg-cell receivers are discussed earlier in this report, 
and characteristics equations are provided for the ideal receiver. The Bragg-cell 
and AO transducer are key producers of first-order distortion effects. Distortion 
effects include AO transducer noniinearities and reflected acoustic signals 
returning from the far end of the Bragg-celL This analysis assumed that acoustic 
signals generated by AO transducer nonlinearities are significantly larger than 
reflected acoustic signals. 

The ideal Bragg-cell receiver contains a single acoustic signal 

=• 2&1 cosjair? ,t +• Eq. 13. 

which results in the diffracted laser information signal given in Equation 12. 
f^is the frequency of the acoustic signal, which is directly related to the 
frequency of the input IF signaL (See Equation 7.) 

AO transducer nonlinearities are producing undesired acoustic signals 
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in the Bragg-cell. Undesired or unwanted acoustic signals are described as 
a cosinusoidal series of N components (harmonics and inter modulation 
components at a significant energy level) which are characterized as: 


N 

u &•)=•!> cos 

ms. Z 

Total unwanted acoustic signals in Bragg-cell. The AO 
transducer will generate harmonics of f^ and other 
intermodulation frequency components. 

Amplitude of a th. unwanted acoustic signal component. 
Frequency of n th. unwanted acoustic signal component. 
Harmonics of f^ and intermodulation frequency components 
are present. 

Phase variations of n th. unwanted acoustic signal component. 
The total acoustic signal in the Bragg-cell is found by summing Equations 13 
and 14: 

wtel* ^ ^ 

Equation 15 can be rewritten, with a lot of work, as: 

— 2 (£) A Ot) cas fair j t + C^O + /S&TJ ***“ 16 ' 


where 

ult')* 


t-W to = 

= 


* 


5 


4- <*„&) 


Eq. 14. 
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where 




t an 


.1 


j> S„(*) 

msg 

i 


K m SlU^A^t +A©^j 

C„ &) a K» COS +• A ©/•) 

q ,(<0 

— Z('b') 


A^ — $«-£ t 
A ©A, =• d i*&')- <*»&) 


The instantaneous frequency variation, resulting from 8 ( t ) phase 
variations in Equation IS, is 
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A?iOfc) = L . 


N 


SJ$ 






AOJmS, 


mat 


— | - M 

*(•<=) + i+Z c ^ 

I Ma ^ 


M 


AW m C m C-t) 


mat 


Eq. 17. 


TT A 2 (-t) 


The instantaneous frequency variation is a very complex function of 
amplitudes, frequencies, and phases of ail unwanted acoustic signal components 
in the Bragg-celL Equation 17 can be rewritten using the Fourier series 
expansion as a fixed frequency plus a series of harmonically related frequency 
components. A first-order approximation is the reduction of Equation 17 to 
a fixed frequency, 

A£ i (-fa') = ■¥■* ** 1 *- 


f 2 is used to approximate the resultant frequency component, which 
is produced by unwanted acoustic signal components propagating in the 
Bragg-cell. 

Resultant phase variations are approximated as: 


y3(t}A a^r^-fc 


Eq. 19. 
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The total acoustic signal (Eq. 16) can be approximated as: 


uj(-t) = -2 It) A(t) c<js£air?,t +air£ J t +«,(*Tj 


Eq. 20. 


The total instantaneous frequency of w( t ) (Equation 20) is: 


?i&) = ?, + -? a 


Eq. 21. 


Equation 12 can be easily changed to approximate unwanted acoustic 
signals. The modified Equation 12 results in the following equation: 

-fi.Cf) = L\(t) cosjair^+^+^t + S'(t)]. ^ 22 . 


The actual displacement of the diffracted laser information signal (Equation 
22) relative to the undeflected laser light is 



V, 


Eq. 23. 


Equation 23 reveals that d is a random variable since it is a summation 
of two random variables f 1 and f 2> Conceptually, the spot size of the 
diffracted laser information signal is increased by the unwanted acoustic signal 
components. The actual displacement d is directly calibrated to a photodiode 
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ceil in the linear detector array, which corresponds to a specified frequency 
bin. 

f 1 and f are assumed to be inaependent random variables for this 
analysis. Pdf characteristics of 1 1 and f are assumed to be either constant 
or Gaussian. The pdf characteristic of d ’ is the convolution of pdf (£ ) 
and pdf ( f ) , if the random variables are independent (Reference 4.) 

U 






pi 4 ($) df, 


Eq. 24. 


where 




Each frequency bin in the linear photodiode array exactly corresponds 
to Ad bin coverage in displacement space. 

The output frequency of the detector array is 


$ 0 = K,d 


Eq. 25. 


Substituting Equation 25 into Equation 23, the result is: 



K a P ($, + S,') 
Ns 


Eq. 26. 


The Bragg-cell receiver's measured output frequency can be stated in 
terms of the Input IF signal's frequency by using Equations 26 and a modified 
version of Equation 7. The output frequency is: 
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Eq. 27. 


■^o- K 3 ($zr + ?a 


where 


K 3 — 


K a F 3- 
K Vs. 


£ 


The equivalent IF frequency of unwanted acoustic signal 

A 


components approximated as acoustic frequency f . 


f is a linear combination of two random variables f __ and 
0 ■ IF 4 , 

, which are assumed to be independent. The pdf (f ) is the convolution 
of pdf (f ) and pdf ( fi* ) , which is expressed mathematically as: 

IF Z 





Eq. 28. 


The pdf (f Q ) characteristic is functionally divided into sampled A£ q 
bins by the photodiode cells in the linear detector array. This frequency domain 
binning process is equivalent to discrete sampling of the pdf ( £ ) 

characteristics in A£ Q bin widths. A£ Q is the basic frequency 
resolution of the Bragg-cell receiver at the linear detector array. The actual 
pdf characteristic of the Bragg-cell receiver's frequency parameter is a discrete 
characteristic, which is a function of the pdf of the IF signal frequency and 
pdf characteristic of the equivalent unwanted signal components. 
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SIGNAL PROCESSING ANALYSIS 


INTRODUCTION 

Signal processing analysis results are presented in this section. The 
Bragg-cell receiver model is implemented on a Hewlett Packard Integral 
computer. Computer software program listings are contained in Appendix 
A of this report. 

Pdf analysis of the Bragg-cell receiver's frequency parameter is 
performed, assuming constant and Gaussian pdf characteristics for the input 
signal's IF frequency. Pdf analysis is performed for selected frequency 
resolutions of the Bragg-cell receiver. Also, statistical analysis are performed 
on output pdf signals and summary results are presented. 

BRAGG-CELL RECEIVER MODEL RESULTS 

The Bragg-cell receiver is math modeled in the previous section. 
First-order model implementation is shown to be a convolution of pdf (f_) 
and pdf (f for the Bragg-cell receiver's output frequency pdf characteristic. 
Equation 28 is the final descriptive equation for the output frequency pdf 
characteristic. 

Bragg-cell receiver signal processing results are summarized in Table 
1. Signal processing analysis programs are presented in Appendix A. Figures 
5 through 70 contain detailed pdf plots for the Bragg-cell receiver's output 
frequency parameters. These plots are made for selected pdf characteristic 
of fjp and f j and selected Bragg-cell receiver frequency resolutions. 

Table 1 defines the input signal or fjp in terms of DFl or DF1 and 
Cj. The pdf (fjp) is assumed to be a constant if only a value for DFl is given. 
DFl (MHz) indicates the frequency excursion of the input frequency (fjp) 

For a Gaussian pdf DFl(MHz) and a^MHz) is selected as 

1.166 MHz or 4.166 MHz for Gaussian pdfs. 
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TABLE 1 


BRAGG-CELL RECEIVER 
SIGNAL PROCESSING RESULTS 
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FIGURE 
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OUTPUT FREQUENCY: 
STATISTICAL RESULTS 

DFI(MHz) 

(MHz) 

DF2(MHz) 

o 2 (MHz) 

o(MHz) 

SKEW 

KURTOSIS 

23 

10 

1.166 

2 


0.5 

1.755 

0.007 

2.850 

24 

10 

1.166 

2 


1.0 

1.755 

0.086 

2.844 

25 

10 

1.166 

2 


2.0 

1.753 

0.111 

2.824 

26 

10 

1.166 

5 


0.5 

2.208 

-0.006 

2.713 

27 

10 

1.166 

5 


1.0 

2.209 

-0.042 

2.714 

28 

10 

1.166 

5 


2.0 

2.207 

0.023 

2.704 

29 

25 

4.166 

15 


1.0 

6.026 

-0.001 

2.624 

30 

25 

4.166 

15 


2.0 

6.026 

-0.008 

2.625 

31 

25 

4.166 

15 


5.0 

6.022 

0.009 

2,617 

32 

25 

4.166 

24 


1.0 

8.124 

0.000 

2.327 

33 

25 

4.166 

24 


2.0 

8.124 

0.003 

2.327 

34 

25 

4.166 

24 


5.0 

8.122 

0.000 

2.324 

35 

10 


2 

0.333 

0.5 

2.915 

0.003 

1.831 

36 

10 


2 

0.333 

1.0 

3.257 

1.174 

2.192 

37 

10 


2 

0.333 

2.0 

2.748 

0.080 

1.799 

38 

10 


5 

0.833 

0.5 

3.012 

-0.004 

1.973 

39 

10 
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0.833 

1.0 

3.102 

-0.019 

1.973 

40 

10 
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0.833 

2.0 

3.055 

0.564 

2.059 
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FIGURE 

INPUT 

SIGNAL 

DISTORTION SIGNAL 

FREQUENCY 

RESOLUTION 

(MHz) 

OUTPUT FREQUENCY: 
STATISTICAL RESULTS 

DFI(MHz) 

mm 

DF2(MHz) 

o 2 (MHz) 

o(MHz) 

SKEW 

KURTOSIS 

41 

10 


9 

1.5 

0.5 

3.255 

0.001 

2.241 

42 

10 
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1.5 

1.0 

3.255 


2.241 

43 

10 


9 

1.5 

2.0 

3.255 


2.239 

44 

25 


5.0 

0.833 

1.0 

7.288 


1.831 

45 

25 


5.0 

0.833 

2.0 

7.306 

-0.165 

1.834 

46 

25 


5.0 

0.833 

5.0 

6.87 


1.799 

47 

25 


15 

2.5 

1.0 

7.653 


2.037 

48 

25 


15 

2.5 

2.0 

7.653 


2.037 

49 

25 


15 

2.5 

5.0 

7.645 

0.063 

2.041 

50 

25 


24 

4.0 

1.0 

8.254 

0.000 

2.281 

51 

25 


24 

4.0 

2.0 

8.254 

0.002 

2.280 

52 

25 


24 

4.0 

5.0 

8.252 

0.006 

2.278 

53 

10 

1.666 

2.0 

0.333 

0.5 

1.683 


2.842 

54 

10 

1.666 

2.0 

0.333 

1.0 

1.683 


2.842 

55 

10 

1.666 

2.0 

0.333 

2.0 

1.680 


2.804 

56 

10 

1.666 

5.0 

0.833 

0.5 

1.845 

-0.003 

2.884 

57 

10 

1.666 

5.0 

0.833 

1.0 

1.845 

-0.003 

2.884 

58 

10 

1 .666 

5.0 

0.833 

2.0 

1.845 

0.002 

2.882 
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The distortion signal or f * is quantified in terms of DF2 or DF2 and 
a The pdf(fb is assumed to be a constant if only a value for DF2 is given. 

U u 

DF2(MHz) indicates the frequency excursion of the distortion signal f 2 * 

For a Gaussian pdf(f ^), DF2(MHz) and o 2 (MHz) indicates six times 
sigma and standard deviation, respectively. DF2 is selected as 2, 5, 9, 15, 
and 24 MHz. a 2 (MHz) is selected as 0.333, 0.833, 1.5, 2.5 and 4.0 MHz 
for Gaussian pdfs. 

Bragg-cell receiver frequency resolution (MHz) is chosen to be 0.5, 
1.0, 2.0 or 5.0 MHz. These frequency resolution selections are readily 
expected with current Bragg-cell receiver technology. 

The last three columns of Table 1 contain statistical results of the 
Bragg-cell receiver's output frequency. Sigma (MHz), skew and kurtosis 
are computed for each output frequency pdf characteristic shown in Figure 
5 through 70. 

The pdftf £ characteristic is a staircased trapezoidal .function for pdf(fjp) 
= constant and pdf(f^) = constant. Figures 5 through 22 show numerous 
examples of pdf(f^) characteristics for selected input signal and distortion 
signal constant pdf characteristics. Pdf(f^) is observed to significantly depart 
from a constant characteristic and approach a staircased triangular 
characteristic as DF2 approaches DF1. 

Figures 23 through 34 reveal output frequency pdf(f^) characteristics 
with input frequency characteristics assumed to be Gaussian pdf and distortion 
signal frequency characteristics assumed to be constant pdf. These pdf 
characteristics are staircased Gaussian characteristics. The pdf characteristic 
departs from the Gaussian characteristics as the Bragg w cell receiver's 
frequency resolution is decreased. 
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Figures 23 through 34 reveal output frequency pdfif^) characteristics 
with input frequency characteristics assumed to be Gaussian pdf and distortion 
signal frequency characteristics assumed to be constant pdf. These pdf 
characteristics are staircased Gaussian characteristics. The pdf characteristic 
departs from the Gaussian characteristics as the Bragg-cell receiver's 
frequency resolution is decreased. 

Figures 35 through 52 reveal output frequency pdf(f^) characteristics 
with input frequency characteristics assumed to be constant pdf and distortion 
signal frequency characteristics assumed to be Gaussian pdf. For most 
example plots the long constant pdf characteristic of the input frequency 
tends to dominate the overall pdf characteristic. The pdf characteristics 
are approximately a staircased trapezoidal characteristic. 

The pdf(f^) characteristic is a staircased Gaussian function for pdf(f^p) 
= Gaussian and pdf(f *) = Gaussian. Figures 53 through 70 show numerous 
examples of pdf(f^) for selected input signal and distortion signal frequency 
characteristics. Pdf(f^) is observed to be close to a Gaussian characteristic 
for all DF2 and DF1 selections. For low frequency resolution in the 
Bragg-cell receiver, the pdf characteristic departs from a Gaussian 
characteristic. 
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CONCLUSIONS 


Bragg-cell receivers can accurately measure frequencies of time 
coincident and time overlapped RF signals, which is a big advantage over 
other receiver types. Bragg-cell receivers are available with bandwidths 
to approximately 1 GHz, and frequency resolutions of 100 KHz to 10 MHz. 

The Bragg-cell receiver's frequency parameter is analyzed and results 
are presented in this report. Probability density function analysis and 
statistical analysis results are presented for this receiver for selected 
frequency resolution capabilities. 

The first order math model is developed for the Bragg-cell receiver. 

The input signal's frequency is assumed to be a constant pdf or Gaussian 
pdf characteristic. The primary distortion signal is generated by nonlinear 
acoustic transducer characteristics. The distortion signal is also assumed 
to have a constant pdf or Gaussian pdf characteristic. 

Staircased pdf characteristics are observed for all pdf frequency 
characteristics. The staircased function is introduced by the frequency 
sampling of the photodiode detector array used in the Bragg-cell receiver. 

The output frequency pdf characteristic is observed to depart from the input 
signal's frequency pdf characteristic for many of the wideband RF signal 
cases and wideband distortion signals produced by the AO transducer. 

Distortion signals with large DF2 values can significantly change the 

output frequency's pdf characteristic. Also, the frequency sampling 

(or detector frequency resolution) at the detector can further distort 

the pdf characteristic. Low frequency resolution at the detector produces 

the most significant distortions. EW emitter classification and EW 

direction finding systems can be affected by Bragg-cell receiver distortions. 

These distortions are especially important for wideband radar emitters. 


99 



Math model validation is recommended for the Bragg-cell receiver. 
Model validations are readily performed using an actual Bragg-cell receiver 
and associated test equipment. 
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APPENDIX A 
BRAGG-CELL RECEIVER 
CONVOLUTION PROGRAM 


10 ! "BRAGG" PROGRAM 

20 ! BRAGG-CELL RECEIVER ANALYSIS 

30 ! CONVOLUTION OF F! ARRAY AND F2 ARRAY 

40 ! THIS VERSION PERMITS STORAGE OF CONVOLVED ARRAY UITH ITS FREQUENCY AXIS 

50 ! STORAGE TAXES PLACE PRIOR TO PLOTTING— A STRING NAMED INFOS IS ALSO 

SO ! STORED— SHOULD CONTAIN MINIMUM INFORMATION QN STORED OATA 

70 ! PLOT BOTH INPUT ARRAYS ON X AXIS 

30 ! PLOT CONVOLUTION ON FREQUENCY AXIS 

SO ! LABEL TOTAL VALUES FOUND IN CONVOLUTION ARRAY 

100 ! ARRAYS MAY BE ANY COMBINATION OF FLAT/GAUSSIAN 

110 ! DECEMBER 133S 

120 ! ISP “DO YOU UANT TO SAVE CONVOLVED ARRAYS?" 3 LINPUT "ENTER Y/N~,XP3 
130 I IF KPSC 5,13 ="Y" THEN 30SUB KEEP 

140 DI3P "ENTER INTEGER LENGTH OF LONG ARRAY" 3 INPUT DF1 
ISO IDISP "ENTER LENGTH OF SHORT ARRAY" 3 INPUT OFZ 
ISO 01 SP "MAKE THESE ENTRIES IN MHZ" 

170 01 SP "ENTER DELTA FREQUENCY FOR THE LONG ARRAY" 3 INPUT DFF1 
130 01 SP "ENTER DELTA FREQUENCY FOR THE SHORT ARRAY" 3 INPUT DFFZ 
130 GF2*DF1/<DFr 1/DFF2) 3 N*CFI*Z 3 SGF 1 *QFF ! * 1000000/ S 3 SGFL0*QFF2* 1 000000/6 
200 01 SP "OFT*" i DF1 1" OFZ*" ; OFZ i" N»"jN 3 NMR1* 1 /< SGF1 *SQR< 2*PI ) ) 

210 OPTION BASE 1 3 RAD 3 NMRLO 3 1 /( SGFLQ*SQR( 2*PI ) ) 3 MU*5O0000000 

220 DIM PX(SOO) ,PH< 500) ,PY< 439 ) .TIC 321 .TISISZJ ,FX(S0O> .XC50O) ,Y(5O0) .INFQSC30] 

230 LINPUT "ENTER *<80 CHARACTERS DESCRIBING DATA", INFO* 

240 REDIM PY< N- I > ,FX< N- 1 ) ,PX( N) ,PH< N) 

250 FOR K-1 TO N ! FILL INPUT ARRAYS UITH ZEROES 
2S0 PX<X)=0 3 PH< K ) *0 


270 NEXT X 

280 OISP “ *****PICX COMBINATION*-****" 

230 DISP " 1 -BOTH ARRAYS FLAT" 

300 OISP “ Z-SHORT FLAT, LONG GAUSSIAN" 

310 OISP " 3-LONG FLAT, SHORT GAUSSIAN" 

320 DISP " 4- BOTH GAUSSIAN" 


330 OISP "ENTER YOUR CHOICE" 3 INPUT PIC 
340 IF PIC=Z THEN SHRT 
350 IF PIC*4 THEN GOSUB ARZ 
360 LSTEP-1 

370 FOR '<*1 TO N\Z ! MAKES FLAT LONG ARRAY 
380 PH<X)*1/<DF1-n 
330 NEXT K 

400 SHRT: I MAKES FLAT SHORT ARRAY 

410 IF PIC=3 THEN GOSUB AR3 
420 SSTEP* 1 

430 FOR X*N/4-DFZ/Z TO N/4+CF2/2 
440 PX< K 1 /( OFZ- 1 ) 

450 NEXT K 

460 IF PIC*Z THEN GOSUB ARZ 
470 CNVLV: 

480 YTOT»0 

430 FOR K*t TO N-J 
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500 Y5UM=0 

510 STL=S00-CFF1 a STPL=500+QFF 1 3 5TZL=< STPL-STL)/< N-1 ) a IVAL=STL+STZL/2 

SZO FOR J*1 TO K 

530 YSUM=YSUM+PX< J)»PHCK-J+1 ) 

540 PY(X)=YSUM/(3TZL* 1000000) 

550 NEXT J 

550 IF K non 50=0 THEN DISP “K=** j K 
570 NEXT K 

530 FOR K= 1 TO N-1 ! MAKE FREQ ARRAY FOR XAXIS 
530 FX< K )=IVAL 
□00 IVAL*IVAL+STZL. 

510 NEXT X 
SZO SUP1Y=0 

530 FOR K= 1 TO N-Z ! INTEGRATE 
540 XS=<FX(X+1 )-FXCX)>* 1000000 
550 YS=< PY( X )+PY< K + 1 ) )/Z 
650 3UMY=XS* YS+SUMY 
570 NEXT X 

580 OISP "SUM 3 " ! SUMY 

590 ? IF KPSt 1,11 =“Y” THEN GOSUB XP1 

700 CHOICE: 1 PICK PLOT 

710 DISP " 1-PLOT 80TH ORIGINAL ARRAYS SAME AXIS” 

7Z0 OISP " Z-PLOT CONVOLUTION ON FREQUENCY AXIS" 

730 DISP " 3-EXIT PROGRAM” 

740 DISP “ 4- STORE DATA" 

750 DISP " ENTER YOUR CHOICE” 0 INPUT CH 

750 IF CH*4 THEN GOSUB XEEP 

770 IF CH#1 THEN 7S0 

780 LINPUT “EXTERNAL PLOTTER? Y/N“.HPS 
730 IF CH=3 THEN BEEP 8 DISP ”OONE” Q STOP 
300 IF CH* 1 THEN GOSUB PLT1 
310 IF CH-Z THEN GOSUB PLT3 

3Z0 PLOTTER IS 1 8 PEN -1 8 GCLEAR 8 LOCATE 30,170,25,30 
330 Y1=0 8 Y2=AMAX( Y)* 1 . 1 8 X1=STL 8 XZ-STPL 
340 START: 

350 SCALE X1.X2.Y1.Y2 

850 FXD Z 8 LAXES < XZ-X1 )/ 10, ( Y2-Y1 >/ 10,X1 , Y1 
370 XAXIS YZ ,< X2-X1 )/10,X1 ,XZ 
830 YAXIS XZ,(Y2-Y1 )/10,Y1 ,YZ 

830 MORDAT: ! ADD MORE DATA SAME SCALE 

900 IF CH=1 THEN 950 
910 IF CH=4 THEN 350 

920 DISP "ENTER BRAGG- RECEIVER FREQUENCY RESOLUTION" 

930 INPUT RES 

940 KRES=INT< RES/< 2*QFF 1/N) ) 

950 FOR K*1 TO TOP 

950 PLOT X( X ) ,Y< X ) 8 IF CH-! THEN t©40 
970 IF CH*4 THEN 1040 
980 FOR L-1 TO XRES 
930 IF K+L>*TCP THEN 1050 
1000 PLOT X< X+L) ,Y( X ) 

1010 NEXT L 
1020 K*X+L-Z 
1030 L* 1 
1040 NEXT K 

1050 IF HPS*“Y” THEN EPLOT 
1050 IF CH* 1 THEN G0SU8 PLT2 
1070 IF CH*4 THEN PL OUT 
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t 080 I MfiGc 

•« 

“CONSTANT 

POF 

FOR 

r 

1030 IMAGE 

“CONSTANT 

POF 

FOR 

r 


OF I , QQ.G. “ 
DF1-" ,00.0," 


CONSTANT POF FOR FZ , CFFZ-“.D0.0 
SAUSSIAN POF FOR FZ . OFF!-" ,DD. D 


1100 IMAGE “GAUSSIAN POF FOR FI, QF!-“,DD.O," CONSTANT POF FOR FZ , GFF2 3 “,Q0. 

0 ," 

1113 IMAGE “GAUSSIAN POF FOR FI. OF 5 «'* , GO. Q, “ SAUSSIAN POF FOR FZ , 0FF2-“,00. 

O," 

1120 T3=“3RAGG-CELL RECEIVER POF ANALYSIS" 8 CSIZE 3.3,1 

1130 LORG S 8 MOVE XZ-< XZ-Xt >/Z , YZ+< YZ-Y1 )/9 9 LABEL USING “K" j T* 

1 140 CSIZE 3 

11S0 MOVE X2-< XZ-X1 )/Z , Y2+< YZ-Y1 )/ 1 3 3 IF PIC-1 THEN LABEL USING 1080 j OFF 1 .OFF 
Z a CSIZE 3 

1160 IF PIC-3 THEN LABEL USING 1090 ! 0FF1 ,DFFZ a CSIZE 4 

1170 IF PIC-Z THEN LABEL USING 1100 5 DFF1 ,OFFZ 8 CSIZE 4 

1130 IF PIC-4 THEN LABEL USING 1110 i DFFl.DFFZ 3 CSIZE 4 


1 1 90 MOVE 
1200 MOVE 
X 1 E3“ 



, Y1-< YZ-Y1 )/S 3 LABEL USING “X“ j “FD IN MHZ" 

, Y2-< YZ-Y1 >/Z 3 DEG 3 LOIR 30 3 LABEL USING “X" 


“Pdf < FD) 


1210 IMAGE “TGTAL-“.CO.DCD 

1 ZZ0 MOVE XI ,Y1-< Y2-Y1 )/S 3 LOIR 0 3 LORG Z 3 LABEL USING 1210 { SUMY 
1230 PLCUT: 

1240 LINPUT “DUMP GRAPHICS? Y/N",DG3 
1230 IF DG3»”Y“ THEN DUMP GRAPHICS 
1230 EPLOT : 


1270 HP3=“N” 


1280 LINPUT “PLOT ON EXTERNAL DEVICE? Y/N",PL£S 

1230 IF PLE$-“Y“ THEN PLOTTER IS 70S 3 PEN 1 3 LOCATE 30,110,23,39 3 SOTO START 
1300 IF CH«4 OR CH-Z THEN CHOICE 
1310 BEE? 150,300 3 01 SP "DONE" 

1320 ENO 


1330 PLT1 : ! PLOT SHORT ARRAY 

1340 REDIM X( N) , Y< N) 

1350 IF CH-1 THEN MULT«t ELSE MULT- 1 00000000 
1 3B0 FOR X-J TO N ! FILL X ANO Y ARRAYS 
1370 XIX )*< 3 Y(X)*PX(X)*MULT 
1380 NEXT X 

1390 XI -1 3 X2-N 3 Y1-0 3 Y2-AMAX1 Y)- 1 . 1 3 TOP-N 3 LINE TYPE 1 
1400 PLOTTER IS 1 8 PEN -1 8 SCLEAR 8 LOCATE 30,170,25,80 
1410 SOTO START 
1420 RETURN 


1430 PLT2: ! PLOT LONG ARRAY SAME SCALE 

1440 T*-“ INPUT ARRAYS" 

1430 IMAGE “CONSTANT ARRAY FOR FI, DF1--.OO.D," CONSTANT ARRAY FOR F2 , 0FF2-" , 

□ 0.00 

14S0 IMAGE "CONSTANT ARRAY FOR FI. 0F1--.D0.0." SAUSSIAN ARRAY FOR F2 , DFF2-" , 
OD.DO 

1470 IMAGE “GAUSSIAN ARRAY FOR FI. OFI-".DD.D," CONSTANT ARRAY FOR F2 , DFF2-" , 
00.00 

1480 IMAGE “SAUSSIAN ARRAY FOR Ft. OFt-“,DO.D,“ SAUSSIAN ARRAY FOR F2 . 0FF2-" . 
00.00 

1490 LORG 5 

1300 MOVE X2-< X2-X1 )/Z , Y2+< Y2-Y1 )/20 8 IF PIC-2 THEN LABEL USING 1470 | 0FF1.DFF 
2 8 CSIZE 4 

1510 MOVE XZ-< XZ-X1 )/2 ,Y2+< YZ-Y1 )/20 8 IF PIC-1 THEN LABEL USING 1450 I DFF1.DFF 
2 8 CSIZE 4 

1520 MOVE X2-( X2-X1 )/2 , Y2+< Y2-Y1 )/Z0 8 IF PIC-3 THEN LABEL USING 14S0 | DFFI.DFF 
2 8 CSIZE 4 
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1530 MOVE XZ-<XZ-X1 >/Z,YZ+< Y2-YD/Z0 a IF PIC=4 THEN LABEL USING 1430 ; DFF1 OFF 
2 a C3IZS 4 

1540 LORG 5 a MOVE X2-<X2-X1 >/2.YZ+< YZ-YD /3 a LABEL USING "X" } TS 
1550 IF CH«! THEN MULT- I ELSE MULT * 1 000000000 

1SS0 MOVE X1-<XZ-X1)/7,YZ-< YZ-YD/2 3 OEG 8 LDIR 30 3 LABEL USING "K“ i "VALUE X 
" ; MULT 

1570 FOR K= 1 TO N 
1580 Y< X )-PH< X ) *MtlLT 
1530 NEXT X 

1S00 CH=4 3 LINE TYPE 4 3 TOP-N 
1S10 SOTO MOROAT 
1SZ0 RETURN 

1 S30 PLT3: ! FILL CONVOLUTION ARRAY 

1640 REDIM X( N- 1 ) , Y( N- ! ) 

1650 MULT= 1 000000000 

1GS0 FOR X»1 TO N-t 

1670 X( K >»FX< X 5 3 Y( K ) a PY( X >*MULT 

1630 NEXT X 

1630 TOP-N-1 

1700 RETURN 

1710 ARZ: ! SUB-GENERA TE SAUSSIAN FOR THE LONG ARRAY 

17Z0 Stt=MU-3*SGF1 a 3tp=MU+3*SGF1 9 STZ-< Stp-Stt )/< N\Z > 8 IVL=Stt 

1730 FOR K»t TO N\2 

1740 PH< X 1-NMR1 *EXPC ( -( IVL-MU ) “2 )/( Z*SSF1 *21) a PH< X )«PH< X >*STZ 
1750 IVL-IVL+STZ 
17B0 NEXT X 

1770 IF PIC-4 THEN AR3 
1780 IF PIC-Z THEN CNVLV 
1790 RETURN 

1800 AR3: ! SU8-SENERATE GAUSSIAN FOR THE SHORT ARRAY 

1810 St-t=MU-3*SGFLO S Stp=MU+3*SGFL0 9 STZ*< Stp-Stt )/( N/4+OFZ/Z-< N/4-0F2/Z ) ) 8 I 
VL*Stt 

1 8Z0 FOR X-N/4-0FZ/2 TO N/4+0FZ/Z 

1830 PX(X)»NMRL0*EXP((r<IVL-MUrZ)/(2*SGFL0*Z)) 9 PX( K )«PX( X )*STZ 
1840 IVL-IVL+STZ 
1850 NEXT X 

I860 IF PIC-3 OR PIC-4 THEN CNVLV 
1870 RETURN 

1880 KEEP: I FILENAME ENTRY 

1890 01 SP “FILENAME MUST BE ENTERED IN UNIX PATH FORM /DI SCNAME/ FILENAME" 

1900 LINPUT “ENTER OESCRIPTIVE FILENAME", F* 

1910 CREATE FS,1 .233*2+89 

19Z0 ASSIGN# 1 TO FS 

1930 PRINT# 1 » N, INFOS, FX< ) ,PY< ) 

1940 ASSIGN# 1 TO • 

1950 SOTO CHOICE 
I960 RETURN 
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APPENDIX A 
BRAGG-CELL RECEIVER 

PDF GENERATION 6 STATISTICAL ANALYSIS PROGRAM 


S 1 8RAGGP0F THIS PROGRAM COMPUTES 3RAGG-CELL RCUR. POF AND STATISTICS 
; 3 ; THIS PROGRAM GETS DATA FROM A C0NU0LU0LUTI0N FILE-QCES CEN MOMS 

23 ! CATA IS THEN PLOTTED WITH LABELLED CENTRAL MOMENTS. INT . ETC 

33 OPTION BASE : 

43 DIM FX<530> , PH (530) .INFQSC80] .TSC753 

53 OISP "PLEASE USE FILENAMES WITH Z DIGITS FOR DF1 * 

S3 LINPUT "FILENAME?" ,FS 
73 ASSIGN* 1 TO -? 

83 READ* 1 ! N 3 N=N- 1 
S3 REDIM FX < N > ,PH<N ) 

' 83 READ* ! f INFOS , FX <>. PH < > 

113 FOR K-i TO N 
1 Z3 FX ( K )— * X 1 K ) + 1 33 
133 NEXT K 
1*3 ASSIGN* 1 TO * 


153 

OISP 3 0 

ISP INFOS 3 013 

P 



183 

FI 3-REVS 

» F3 ) 3 

PS-* 1 SC 1 

,1 1 



173 

IF PS- "A 

* THEN 

PIC-1 




180 

IF ps— "3 

" THEN 

PIC-2 




1 80 

IF PS-“C 

“ THEN 

PIC-3 




233 

IF PS-*D 

“ THEN 

PIC-4 




*■* 1 /** 
L, i V 

DFi=VAL< 

F3C Z ,4 

j / « urUu 

-VALCFSC 

s 

V 

, / i / 


OISP "ENTER RE 

SOLUTION 

IN MHZ" 

3 

INPUT 


XST-INTI 

RES/ (2 

•QF1/N) ) 





243 CFM*AMAX< FX 5 -AM INC FX 5 9 FSTPM-QFM/C N-f >*XST 9 OF-FSTPM* t 330000 9 SUMZ-0 
253 3UMY=3 9 3UMS-3 9 SUM3-0 9 SUM4-0 
2S3 FOR K- 1 TO N-l 

273 XS-(FX'K+1 >-FX< K >>* 1 330000 9 YS-<PH<K >+PH<K+1 > )/2 
283 3UMZ=SUMZ+XS*YS 
283 NEXT K 

303 FOR K=! TO N STEP XST 

313 FX<K >=FX(K >*1000030 

323 PH<K)-PH(K >/SUMZ 

333 SUMY=SUMY+PH< K )*FX.< K >*OF 

340 NEXT K 

353 XMEAN-SUMY 

353 FOR K-1 TO N STEP XST 

370 3UMS-SUMS+< FX< K 5-XMEAN )"2*PH( K )*QF 

380 NEXT K 

380 UNC-SUMS 9 3ISMA«SOR< UNO > 

400 FOR K-1 TO N STEP XST 
410 SUM3*SUM3+< FX ( K 5-XMEAN )"3*PH< K )*DF 
423 SUM4=SUM4+< FX< K 5-XMEAN )"4*PH< K )*OF 
433 NEXT K 

440 CM3-SUM3 9 CM4-SUM4 

453 SKEW-CM3/VNC" 1 .5 9 KURT-CM4/UNC"2 

4S3 PLOTTER IS 1 9 PEN -1 9 GCLEAR 9 LOCATE 30,173,25,30 

470 Y1 *0 9 Y2-AMAX1PH 5*100000000*1 .1 9 X1-FX1 1 5/1000000 9 X2»XS*N/ 1 000000/2+530 
480 START: 
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= G3 SCALE XI ,X2 ,Y1 ,Y2 

=33 FXQ 3,2 3 '-AXES ! X2-X 1 >/ 1 0 , < Y2-Y1 )/ 1 3 ,X 5 , Y 1 

'\nAi: l„ , \ A4 *a I / / IU,AI ,A4 

=23 VAX IS X2 Y2-V1 )/18,Y1 ,Y2 

=33 MGRDAT : ! AQQ MORE CATA SAME SCALE 

=43 FOR K= 1 TO N STEF X3T 



IF K+XST^ N THEM 

C2/* 

:ou 









:cv 

rLu 1 r a \ rv > / i wo wo*} ,r 

H< K 

)* t 

4% AAAAAAA •% r* t 

y rLu i 

F\# / M 1 \# 4» v / 

r a^taj i // 

4 (Vt/kA 

1 

03 ,?H< K 

) * 1 83303833 

373 

=23 

NEXT K 
IMAGE 

"CONSTANT 

PDF 

FOR 

Ft 

, OF! *“ 

, ODD . D , * 

CONSTANT 

PDF F 

CR F2 , 

0FF2-" .DDD.O 

D 

= 33 

IMAGE 

•• A»«*f 

i ni* t 

ruF 

run 

F 1 

, DF1 =* 

« 4 

,JUU . u T 

GAUSSIAN 

PDF F 

uH r^ , 

DFF2=" , DDD.O 

U 

*%*% 

a'O'u 

rn 

IMAGE 

“GAUSSIAN 

PDF 

FOR 

FI 

. OF 1 =" 

rs « 

,uuu.u f 

CONSTANT 

PDF 

FOR F2. 

0FF2-" ,DDD . 

uu 

313 

IMAGE 

"GAUSSIAN 

POF 

FOR 

FI 

, DF1 *" 

.ODD.O , " 

GAUSSIAN 

PDF 

FOR F2, 

CFF2-" ,DDD . 


sza 

333 

€43 


73= “ 3R AGG-CELL RECEIVER PDF ANALYSIS" 3 CSIZE 3 
LORE = 3 MOVE XZ-< X2-X 1 5/2 ,YZ+< Y2-YD/9 9 LABEL 


• 3,1 
USING 


v:ua u 


3E3 MOVE X2-< X2-X 1 )/2 ,Y2+< Y2-Y1 >/13 3 IF PIC=1 THEN LABEL 

VOi ** 


333 IF RIC=3 THEN LABEL USING ESQ i DF1 ,DFL0 3 CSIZE 4 

373 IF PIC=2 THEN LABEL USING 333 I 0F1 ,DFL0 3 CSIZE 4 

333 IF PIC=4 THEN LABEL USING 513 s CF1 ,0FL0 9 CSIZE 4 

333 MOVE X2-1X2-X1 )/2 ,Y1 -( Y2-Y1 >/S 3 LORG S 3 LABEL USING 

733 LORG S 3 CSIZE 4 3 LORG 5 


“X" i 
USING 


“K" » 


353 ( CF1 ,DFL0 3 


"Fd IN MHZ" 


713 MOVE X1-<X2-Xn/7.Y2-<Y2-Yl )/2 9 DE3 3 LDIR 93 3 LABEL USINS "K" j "Pdf(Fd) 
X 1 E3" 

723 LDIR 3 3 RAD 
733 IMAGE "INTEGRAL 3 " ,Q.DOD 
743 IMAGE "XMEAN=" .DDDDD.ODO 
7S3 IMAGE "SIGMA*" .ODD. ODD 
7B3 IMAGE "SKEV-" .SDD.OOD 
773 IMAGE "KURT03IS-" .DD.DDO 
733 CSIZE 2.5 3 LORG 2 

7S3 MOVE X1+(X2-X1 )/25 ,Y2-< Y2-Y1 >/10 9 LABEL USINS 733 j SUMZ 

833 LABEL USING 743 j XMEAN/ 1333033 

810 LABEL USING 7S3 » SIGMA/ 1 000000 

823 LABEL USINS 7S0 j SKEW 

333 LABEL USINS 770 ? KURT 

843 LABEL USINS "K" » "RESOLUTION-" »RE3 i * MHZ" 

853 LINPUT "DUMP GRAPHICS? Y/N* ,SDMS 
363 IF SOMSI t ,1 >"Y" THEN DUMP GRAPHICS 
873 IF EPL3=*Y" THEN 900 

330 LINPUT "PLOT ON EXTERNAL DEVICE? Y/N" .ERLS 

393 IF EPL*=“Y" THEN PLOTTER IS 70S 9 PEN t 9 LOCATE 30.113,29.39 9 SOTO START 
933 SEEP 383,403 9 OISP "DONE* 

913 ENO 
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